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ABSTRACT 

o 

The  theories  of  Mosby  and  Amot  on  the  response  of  bound¬ 
ary  layer  properties  are  applied  to  the  computation  of  air 
temperature  and  water  vapor  pressure  at  a  height  of  10  m. 
above  the  sea  surface.  The  effects  of  different  empirical 
constants  on  the  equilibrium  differences  of  the  sea  and  air 
temperatures  and  water  vapor  pressure,  were  investigated 
numerically.  The  computations  indicate  that,  for  all  prac¬ 
tical  purposes,  the  equilibrium  of  sea  and  air  temperatures 
is  reached  in  about  five  hours  when  the  rate  of  change  of 
sea  surface  temperature  along  the  trajectory  of  the  air  re¬ 
mains  constant.  The  numerical  value  of  equilibrium  differ¬ 
ence  depends  on  the  rate  of  change  of  sea  surface  tempera¬ 
ture  along  the  trajectory  of  the  air.  The  results  of  the 
computations  are  presented  in  graphs.  The  application  of 

O 

the  Amot-Mosby  theory  provides  a  method  for  synoptic  numeri¬ 
cal  analysis/forecasting  of  surface  air  temperature  and 
water  vapor  pressure  over  the  oceans .  The  values  of  the 
empirical  constants  should,  however,  be  determined  more 
accurately  from  a  properly  devised  measurement  program  in  a 
variety  of  conditions . 

In  the  numerical  simulation  of  the  feedback  effects, 
the  surface  wind  was  allowed  to  blow  across  an  oceanic 
front,  both  from  warm  to  cold  and  from  cold  to  warm  sides. 

The  surface  air  temperature,  water  vapor  pressure  of  the 


Vll 


air,  and  various  heat  exchange  components  were  computed. 

The  resulting  changes  of  the  properties  of  surface  air  mov¬ 
ing  from  the  cold  to  the  warm  side  are  essentially  quite  dif¬ 
ferent  than  in  air  moving  in  the  opposite  direction.  There 
is  a  characteristic  "overshooting"  of  air  temperature  and 
water  vapor  pressure  at  the  edges  of  the  oceanic  front.  The 
overshooting  is  physically  and  numerically  explainable  on 
the  basis  of  the  dependence  of  equilibrium  differences  (sea 
surface  and  air  temperature)  on  the  rate  of  change  of  sea 
surface  temperature  along  the  trajectory  of  the  air. 

The  computed  air  temperature  and  heat  exchange  compo¬ 
nents  across  the  oceanic  front  are  used  in  the  Sawyer-Bushby 
equation  and  the  thermodynamic  equation.  The  effects  of 
feedback  (at  the  oceanic  fronts)  on  the  state  and  motions  in 
the  lower  layers  of  the  atmosphere  by  different  types  of  air 
flow  were  numerically  simulated.  These  numerical  experiments 
explain  several  weather  and  oceanographic  phenomena  observed 
at  the  fronts .  It  can  be  concluded  that  the  oceanic  fronts 
are  a  significant  feature  in  the  determination  of  the  mean 
position  of  the  atmospheric  fronts  in  the  proximity  of  oceanic 
polar  and  artic  fronts . 


% 
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1.  PURPOSES  OF  THE  STUDY 


The  exchange  processes  between  the  sea  and  the  atmo¬ 
sphere  are  the  primary  cause  for  modification  of  surface 
layers  of  air  over  the  ocean.  In  any  energy  exchange  pro¬ 
cess,  the  differences  between  the  properties  of  the  surface 
layer  of  the  sea  and  the  corresponding  properties  of  the 
surface  layer  of  air  determine  greatly  the  direction  and 
intensity  of  the  exchange  processes. 

Thus,  the  first  purpose  of  this  paper  is  to  study  nu¬ 
merically  the  rate  of  modification  of  the  properties  of 
surface  air  moving  over  the  ocean  and  to  establish  and  test 
numerical  methods  for  synoptic  computation  of  surface  air 
temperature  and  water  vapor  pressure. 

The  horizontal  gradients  of  heating  and  cooling  of  sur¬ 
face  air  masses  create  pressure  gradients  which  affect  the 
movement  of  air.  These  horizontal  gradients  of  heating  and 
cooling  are  expected  to  be  most  pronounced  where  the  under¬ 
lying  surface  has  sharp  thermal  gradients  (i.e.  at  oceanic 
fronts) .  Surface  winds  in  excess  of  twice  the  gradient 
wind  are  often  observed  at  oceanic  frontal  regions  (Cormier 
and  Kindle,  1971).  It  is  of  further  interest  to  point  out 
that  the  mean  position  of  atmospheric  polar  and  arctic  fronts 
over  the  oceans  is  in  the  vicinity  of  the  oceanic  polar  and 
arctic  fronts.  Consequently,  the  second  objective  of  this 
paper  is  to  report  on  a  preliminary  numerical  study  of  the 


effects  of  oceanic  (polar)  fronts  on  the  atmosphere  and  to 
study  the  effects  of  exchange  processes  on  other  frontal 
phenomena  such  as  mesoscale  local  winds  and  fogs . 

Use  will  be  made  only  of  properties  which  are  measured 
during  routine  (synoptic)  weather  observations  or  which  can 
be  easily  derived  from  these  observations  so  that  the  re¬ 
sults  can  be  directly  applied  to  routine  weather  and  oceano¬ 
graphic  analyses/forecasts. 


2.  THEORY  OF  THE  LARGE-SCALE  RESPONSE  OF  SURFACE  AIR 
TEMPERATURE  AND  WATER  VAPOR  PRESSURE  TO  THE 
CORRESPONDING  PROPERTIES  OF  THE  SEA  SURFACE 

The  various  theories  of  the  response  of  maritime  air  to 
the  properties  of  the  sea  surface  were  reviewed  by  Roll 

O 

(1965) .  A  Norwegian  theory  originally  reported  by  Amot  (1944) 
and  later  by  Mosby  (1957)  and  B^yum  (1962)  is  presented  here 
with  a  few  modifications.  This  theory  is  adaptable  to 
numerical  computations  using  parameters  available  in  routine 
meteorological  observations.  In  this  theory  the  trajectory 
of  the  surface  air  is  followed.  The  surface  air  properties 
generally  refer  to  a  standard  height  of  10  m.  above  the  sea 
surface.  The  detailed  profile  of  these  properties  with 
height  (lapse  rate)  at  any  given  time  is  not  specified  as 
it  is  not  available  from  routine  synoptic  observations .  The 
theory  can  be  considered  as  one  of  turbulent  transfer  in 
which  the  intensity  of  turbulent  diffusion  is  related  to 
the  "transfer  potential"  (the  difference  between  the  prop¬ 
erties  of  the  sea  surface  and  those  of  the  air)  and  wind 
speed. 

An  air  parcel  moving  over  the  sea  surface  has  its  tem¬ 
perature  changed  mainly  by  the  convective  turbulent  trans¬ 
fer  of  sensible  heat.  If  the  temperature  of  the  sea  surface 
(T  )  changes  along  the  trajectory  of  the  air,  the  corre¬ 
sponding  individual  change  of  sea-air  temperature  difference 
AT  =  Tw  -  Ta  is : 


dAT  dTw  dTa 


dt  dt  dt 


(1) 


where  Ta  is  the  air  temperature. 

The  change  due  to  sea-surface  temperature  effects  along 
the  trajectory  is  essentially  composed  of  two  terms: 
dTw  9Tw 


dt 


+  \V  •  V  T 


at 


w 


(2) 


The  "local  change"  of  sea  surface  temperature 


9Ta 
9 1 


is  small 


(usually  much  smaller  than  0.01°  C./hour)  and  can  be  neglect¬ 
ed  especially  in  comparison  with  the  last  term  (  \V  •  VTW) , 
which  represents  the  sea  surface  temperature  change  along 
the  trajectory  of  the  surface  air.  \V  is  the  velocity  of 
the  surface  air  parcel  and  can  be  assumed  to  be  the  surface 
wind  speed.  Thus  equation  (1)  reduces  to: 
dAT  dTa 


dt 


=  W  •  VTW  - 


dt 


(3) 


Synoptic  sea  surface  temperature  analyses  as  well  as  analyses 
of  surface  wind  are  available  routinely. 

dT 

The  individual  change  of  air  temperature  ^-a.  can  be  ob¬ 
tained  from  the  first  law  of  thermodynamics: 
dT  dQ  R  T  dp 

_ 1  =  _  +  a  (4) 

dt  cpdt  Cp  p  dt 

where  R  is  the  gas  constant,  p  the  pressure,  Cp  the  specific 
heat  of  the  air,  and  Q  the  addition  of  heat.  The  last  term 


4 


includes  adiabatic  cooling  of  air  by  advection,  the  magni¬ 
tude  of  which  can  be  computed  for  various  conditions .  This 
computation  has  been  performed  by  Amot  who  found  that  ad- 
vective  adiabatic  cooling  of  the  surface  air  is  insignifi- 

#  O 

cant  for  moderate  wind  velocities.  At  55°  N.  Amot  found 

that  if  IS-  is  as  large  as  1  mb .  per  hour,  then 
3t 

RTa9p 

- =  0.078°  C  per  hour. 

cpP9t 

In  the  present  sea-air  interaction  problem  the  direct  radia¬ 
tive  heating  and  cooling  of  the  air  will  be  neglected. 

The  convective  transfer  of  sensible  heat  is  a  function 
of  AT  (see  e.g.  Laevastu,  1960) .  Furthermore,  Mosby  (1933) 
and  B0yum  (1962)  found  from  observations  that  the  changes  of 
air  temperature  can  be  well  represented  as: 

dTa 

-  =  K  (Tw  -  Ta)  -C  (5) 

dt 

where  K  and  C  are  positive  constants.  Substituting  (5)  into 
(3),  re-arranging,  and  assuming  \V  •  VTW  to  be  constant, 
equation  (3)  can  be  transformed  into  the  integrable  form: 


dAT 

- +  K  AT  =  C  +  \V  •  VT  ( 6 ) 

dt 


5 


which  after  integration  gives : 


AT 


AT  e 

O 


-Kt 


dx  ) 


(7) 


when  r  is  distance.  This  is  a  useful  formula  for  computation 
of  AT  where  ATo  is  at  the  time  Tq .  The  empirical  constants 
of  C  and  K  must  be  determined  from  observational  data. 

We  can  also  derive  a  useful  formulation  for  numerical 
computation  of  surface  air  temperature  over  the  ocean  if  we 
assume  as  before  that  the  temperature  of  an  individual  parcel 
of  air  is  changed  proportionally  to  the  sea-air  temperature 
difference  (an  assumption  analogous  to  Mosby's  (1933)  empiri¬ 
cal  results) : 

dT 

- 5-  =  K  (Tw  -  Ta)  -  C  (5a) 

dt 


Expanding  the  term  on  the  left,  the  equation  becomes: 


8T. 


dt 


+  \V  •  ATa  =  K  (Tw  -  Ta)  -  C 


(8) 


Equations  (8)  and  (5)  are  equivalent.  Equation  (6)  and  (7) 

are  frequently  used  for  research  purposes  and  equation  (8) 

is  normally  applied  operationally. 

Equation  (8)  can  be  expanded  to  another  form.  Assume 

that  the  local  sea  surface  temperature  does  not  change 

with  time:  — =  0  and  add  this  zero  term  to  the  left  side. 
9 1 


6 


Add  the  term  W  •  VTW  to  both  sides  of  the  equation: 


9T, 


w 


9T 


a 


9t  9t 


+  \V  •  VTW  -  W  •  VTa  =  “K(TW  -  Ta)  +  C 


+  W  •  VTW 


(9) 


+  W  •  V(TW  -  Ta)  +  K  (Tw  -  Ta)  -  C 


9t 


W  •  VT, 


w 


(10) 


d<Tw  -  Ta> 


+  K(TW  -  Ta)  -  C  =  W  •  VTW 


(ID 


dt 


The  same  argument  as  used  for  air  temperature  change 
can  be  used  for  water  vapor  pressure  (ea)  change  over  water. 
The  theory  might  be  somewhat  more  valid  for  ea  as  there  are 
no  appreciable  radiation  effects  on  ea.  The  empirical  con¬ 
stants  (C  and  K)  must,  however,  be  different  for  ea.  The 
values  of  the  constants  reported  in  the  referenced  literature 
are  summarized  in  Table  1.  Obviously,  additional  investiga¬ 
tions  are  desirable  for  determination  of  accurate  values  for 
these  constants. 

Other  theories  and  approaches  for  prediction  of  air 
temperature  changes  over  the  water  were  also  evaluated. 
Considering  the  heat  content  of  a  column  of  air  moving 
across  the  sea  surface,  Burke  (1945)  derived  an  integral 
formula  (which  he  evaluated  graphically  in  a  series  of 
charts)  to  predict  the  temperature  of  the  air  as  a  function 
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Values  of  Constants  K  and  C  in  Equation  (5)  as  Reported 
in  the  Referenced  Literature. 
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of  initial  air  temperature,  length  of  trajectory  over  the 
sea,  initial  lapse  rate,  and  average  sea  surface  tempera¬ 
ture.  He  then  used  the  charts  to  predict  air  temperature 
in  the  transformation  of  polar  continental  air  to  polar 
maritime  air.  Although  the  predictions  close  to  the 
eastern  coasts  of  the  continents  (where  large  changes  are 
taking  place)  were  relatively  good,  his  charts  were  rather 
difficult  to  use  over  the  major  parts  of  the  open  oceans  — 
especially  at  lower  latitudes.  The  charts  also  underesti¬ 
mated  the  final  air  temperature  in  these  areas .  Further¬ 
more  the  formulas  are  not  well  suited  for  operational  com¬ 
puter  work,  primarily  because  of  difficulties  in  obtaining 
the  correct  initial  lapse  rate  near  the  surface  and  in 
determining  a  realistic  trajectory  of  the  air. 

Another  formula  for  estimating  the  air  temperature 
from  the  surface  temperature  has  been  proposed  by  Duquet 
(1960) .  The  formula,  although  correct  and  elegant,  is  quite 
cumbersome  to  use.  Its  simplified  form  requires  the  use  of 
a  lapse  rate  so  that  his  formula  leads  to  the  same  diffi¬ 
culties  as  encountered  by  Burke. 


3.  THE  SPEED  OF  RESPONSE  TO  Ta  AND  ea  TO  T„  AND  etI 

a  a  w  w 

AND  THE  AIR-SEA  EQUILIBRIUM  VALUE 

Equation  (7)  is  used  for  computation  of  the  change  of 
AT  with  time  when  air  is  moving  over  the  water  using  various 
combinations  of  C  and  K  and  various  rates  of  sea  surface 

temperature  change  (Vr  iZw)  designated  hereafter  as  ATW. 

8r 

Some  of  the  results  are  shown  in  Figures  1  through  4. 

The  difference  AT  =  Tw  -  Ta  (and  corresponding  ew  -  ea) 
reaches  an  "equilibrium  value"  (for  most  practical  purposes) 
after  about  5  hours  and  this  equilibrium  value  is  indepen¬ 
dent  of  the  initial  AT ,  but  dependent  on  ATW  and  on  the  con¬ 
stants  C  and  K  (compare  Figures  1  and  3) .  The  equilibrium 
values  with  different  ATW  and  different  constants  of  C  and 
K  are  shown  in  Figure  4.  The  observed  equilibrium  values 
of  AT  and  Ae  at  Weather  Ship  M  (from  B0yum,  1962)  with 
different  ATW  and  Aew  (i.e.  in  this  case  with  different 
wind  directions  and  speeds)  are  shown  for  comparison  in 
Figures  5  and  6  respectively. 

As  there  is  still  some  uncertainty  about  the  most  accu¬ 
rate  C  and  K  values  to  be  used  in  synoptic  numerical  com¬ 
putations,  it  was  found  desirable  to  demonstrate  the  effects 
of  changing  C  with  constant  K  as  well  as  changing  K  with 
constant  values  of  C.  The  results  of  this  computation  are 
given  in  Figure  7. 
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Figure  2.  Change  of  AT  With  Time  and  With  Different  AT, 

Initial  AT=2.5°  C. 
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Figure  3.  Change  of  AT  With  Time  and  With  Different  AT, 

Initial  AT=1.25°  C. 


CNJ 

CD 


^  & 

“  ZD 

o  hr 


Q.  < 


IS 

GO 

O  ll» 
CD 

LU 

CD 


CD 


?  £ 


CNJ 


CD 

I 


cc. 

ZD 

CD 

» — I 

QZ 

LU 

Q- 

CD 

o 


H3AVM1  SyriOH  31  M3idtf 

3o  vi  -  Hi  wniyamnoH 


HD 

c 

a 

4-J 

o 

0) 

-p 

c 

aJ 

-P 

cn 

c 

o 

u 

-P 

a 

<u 

p 

a; 

4-J 


jG 

«P 


E-» 

<3 


P 

D 


Eh 

< 

-P 

C 

<U 

P 

a) 

<p 

<P 

*H 

Q 

JG 

4-> 


•H  5 


nd 

C 


O1  m 

W 

4H 

O 

C 

O 

w 

*H 

P 

<0 

cu 

s 

o 

u 


0) 

p 

3 

Cn 

•H 


15 


.  vi  -  Mi  wniNaiunoa 


ISO, 

9 


CSJ 

O 


O 


CSJ 


9 


NS 

I 


16 


RATE  OF  CHANGE  OF  SEA  TEMPERATURE 
°C  PER  1  HOUR  TRAVEL  (aT„) 

Figure  5.  Equilibrium  Values  of  AT  With  Different  AT  Observed  by 

B0yum  (1962) .  W 


aw  v3  -  Ma  wniyaninoa 
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RATE  OF  CHANGE  OF  SURFACE  VAPOR  PRESSURE 
mb  PER  1  HOUR  TRAVEL  (aew) 

Figure  6.  Equilibrium  Values  of  Ae  With  Different  Aew  Observed  by 

Bjrfyum  (1962) . 
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Figure  7.  The  Effect  of  Changing  C  With  Constant  K  on  the  Equilibrium 
Values  and  the  Effect  of  Changing  K  With  Constant  C  on  the 

Equilibrium  Values. 
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Figures  1  through  3  indicate  that  when  AT  (i.e.  Tw  -  Ta) 
reaches  an  equilibrium  value,  the  turbulent  transfer  co¬ 
efficient,  which  is  proportional  to  AT,  becomes  a  function 
of  ATw.  This  is  evident  in  the  increase  of  equilibrium  AT 
with  increasing  values  of  ATW. 

It  is  of  interest  to  point  out  that  when  air  moves  to¬ 
wards  colder  water  (negative  ATW) ,  the  difference  between 
the  sea  and  air  temperature  (AT)  decreases.  This  has  been 
observed  in  nature  and  can  be  explained  physically  by  the 
increase  of  stability  of  the  lower  layers  of  the  air  and 
consequent  supression  of  turbulent  transfer  caused  by  cooling 
of  the  air  close  to  the  surface.  The  evaporation  rate  is 
suppressed  in  the  same  manner.  There  must  also  be  a  time 
lag  in  the  adjustment  of  the  properties  of  the  surface  air 
to  the  properties  of  the  sea  surface. 

When  the  air  moves  towards  warmer  water  (positive  ATW) 
then  AT  increases.  This  is  physically  explainable  by  the 
increase  of  turbulence  caused  by  heating  from  below. 

Numerical  synoptic  analysis/forecasting  of  Ta  and  ea 
has  been  carried  out  at  Fleet  Numerical  Weather  Central, 
Monterey  from  1965  to  1969  using  the  method  of  equation  (8) 
(see  also  Carstensen  and  Laevastu,  1966) .  Some  verifi¬ 
cations  of  these  analyses  are  shown  in  Figures  8  and  9. 

Amot  (1944)  presented  the  argument  that  the  heat  ex¬ 
change  between  the  sea  and  the  air  should  be  zero  at  the 
equilibrium  value  of  AT  when  ATW  equals  0  (i.e.  the  air 
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Figure  8.  Difference  Between  Reported  Air  Temperature 
Minus  Analyzed  Air  Temperature  -on  15Z,  18Z  and 
and  21Z,  28  March  and  00Z  29  March  1966. 
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Figure  9.  Difference  Between  Reported  Water  Vapor 
Pressure  of  the  Air  Minus  Analyzed  Water  Vapor 
Pressure  on  15Z,  18Z  and  21Z,  28  March 
and  00Z  29  March  1966. 
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moves  along  sea  surface  isotherms) .  Using  Mosby's  measure¬ 
ments  (1933  and  1942)  at  two  different  heights  above  the  sea 
(6  and  24  m. ) ,  Amot  found  that  the  difference  in  temperature 
was  0.01°  C.  per  meter,  which  equals  the  dry  adiabitic  lapse 
rate.  The  heat  exchange  is  assumed  to  be  zero  under  dry 
adiabatic  conditions  if  the  vertical  transport  of  heat  by 
direct  conduction  and  radiation  is  small  compared  to  turbu¬ 
lent  transport.  If  t  increases,  then  equation  (7)  tends  to 


C  1 

become  AT  =  —  +  —  V. 


3T, 


w 


therms 


3T, 


K  K 


3r 


If  air  moves  along  the  iso- 


3r 


—  =  0  or  it  is  calm  (Vr  =  0) ,  then  the  term  con¬ 


taining  these  values  becomes  0  and  AT  =  —  in  equilibrium 

conditions.  In  these  circumstances,  equation  (5)  takes  the 


dT 


a  


form 

dt 

the  air  occurs. 


=  0  and  no  heating  or  cooling  of  the  lower  layer  of 


Mosby  measured  AT  data  at  heights  of  6  and  24  meters  and 

3Tt 


found  that  AT6=  0.46+  3.6  V 

3T, 


V, 


Lw 


„  -H-  and  AT24  =  0.64  +  8.9 
r  3r  24 

3Tt 


3r 


respectively.  Assuming  V  — JsL  to  be  0  under  equilib- 

r  3r 

o 

rium  conditions  as  mentioned  above,  Amot  concluded  from 

Mosby's  measurements  that  heat  exchange  between  the  sea  and 

C 

atmosphere  should  be  zero  when  AT  =  —  =  0.40  +  O.Olh, 

r 

where  h  is  in  meters .  This  assumption  postulates  that 


K 


changes  linearily  with  height. 

B0yum  reasoned  that  the  heat  exchange  may  be  zero  and 
equilibrium  conditions  may  prevail  even  when  a  temperature 
difference  in  excess  of  the  adiabatic  rate  exists  between 
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the  sea  surface  and  the  air.  The  sea  surface  temperature  is 
normally  measured  by  bucket  or  intake  thermometers  and  these 
methods  do  not  measure  the  real  temperature  of  the  cooler 
surface  film.  A  cooler  surface  film  is  known  to  exist  in 
laboratory  conditions  and  is  caused  by  evaporation  and  sen¬ 
sible  heat  loss  from  the  surface  film.  Furthermore,  the 
sea  surface  temperatures  measured  with  intake  thermometers 
are  known  to  have  a  positive  bias. 

The  problem  of  the  cooler  surface  film  has  not  been 
fully  explored  in  natural  conditions  and  the  film  tempera¬ 
ture  in  various  wind  and  sea  surface  roughness  is  not  known 
with  any  reasonable  accuracy. 

It  should  also  be  emphasized  that  an  equilibrium  con¬ 
dition  does  not  necessarily  mean  that  no  exchange  takes 
place  between  the  sea  and  the  atmosphere.  The  loss  of  heat 
from  the  sea  can  equal  the  removal  of  heat  from  the  atmo¬ 
sphere. 
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4.  EFFECTS  OF  AIR  FLOW  ACROSS  OCEANIC  FRONTS 


In  a  numerical  study  (to  be  published)  of  the  effects 
of  sea  surface  temperature  anomalies  on  the  feedback  of 
energy  from  the  sea  to  the  atmosphere,  it  became  obvious 
that  the  greatest  effects  on  the  atmosphere  must  be  found 
at  oceanic  fronts.  (For  a  summary  on  the  oceanic  fronts 
see  Laevastu  and  La  Fond,  1970.) 

In  the  numerical  simulation  of  the  oceanic  frontal 
effects  using  equation  (7) ,  an  oceanic  front  of  4°  C. 
breadth  (16°  to  20°  C.) ,  a  0.5°  C./10  km.  gradient,  and 
a  10  km.  computational  grid  size  were  selected.  A  surface 
wind  of  5.55  m./sec.  was  programmed  to  blow  from  the  cold 
to  the  warm  side  in  one  experiment  and  vice  versa  in  a 
second.  The  saturation  water  vapor  pressure  of  the  sea 
surface  was  computed  from  the  prescribed  sea  surface  tem¬ 
perature  (Tw)  and  the  surface  air  temperature  (Ta) .  Water 
vapor  pressure  (ea)  was  computed  using  equation  (7)  and 
coefficients  of  C  =.15,  K  =.21  for  Ta  and  C  =.18,  K  =  15 
for  ea.  Some  of  the  results  are  given  in  Figures  10  through 
13. 

When  the  wind  direction  was  from  the  warm  toward  the 
colder  side  (Figures  10  and  11)  the  differences  of  AT  = 

(Tw  -  Ta)  and  Ae  =  (ew  -  ea)  decreased  at  the  southern  edge 
of  the  oceanic  front  and  increased  at  the  northern  edge. 

Thus  an  apparent  "overshooting"  of  At  and  Ae  occur  at  the 
edges  of  the  front. 
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Figure  11.  Change  of  ea  Across  a  Front  After  Wind  has  been 

Blowing  From  Warm  to  Cold  Side  1  Hour. 


When  the  air  moves  from  the  cold  to  the  warmer  side 


(Figures  12  and  13)  AT  and  Ae  increase  at  the  north  side  of 
the  front  and  decrease  at  the  southern  edge. 

"Overshooting"  at  the  edge  of  the  sea  surface  tempera¬ 
ture  gradient  is  actually  an  adjustment  process  of  equilib¬ 
rium,  AT  to  various  rates  of  changes  of  ATW-  along  the  tra¬ 
jectory  (see  Figures  1  through  3) .  As  the  air  on  the  warm 
side  in  the  area  of  weak  temperature  gradient  moves  across 
the  front,  ATW  becomes  negative  and  AT  decreases.  At  the 
cold  side,  ATW  becomes  weaker  and  AT  increases.  The  "over¬ 
shooting"  in  air  moving  from  the  cold  to  the  warm  side  can 
be  similarly  explained.  Several  consequences  of  this  "over¬ 
shooting"  result,  some  of  which  are  explored  in  Section  5. 

The  response  of  the  surface  air  to  the  properties  of 
the  sea  surface  with  respect  to  the  direction  of  air  flow 
across  the  oceanic  front  cause  differential  heating  and 
cooling  of  the  air  and  affect  the  heat  exchange.  The  com¬ 
putation  of  sensible  and  latent  heat  was  accomplished  with 
the  following  formulae  (Laevastu  1960) : 

Qh  =  39  (0.26  +  0.077V)  (Tw  -  Ta)  (13) 

Qe  =  (0.26  +  0.077V)  (ew  -  ea)  Lt  (14) 

where  is  sensible  heat  exchange  (g  cal.  cm.-2  (24  h)-1), 
Qe  is  latent  heat  exchange,  V  is  wind  speed  (m.  sec.-1)/ 
and  L.J.  is  the  heat  of  vaporization. 
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Figure  12.  Change  of  Ta  Across  a  Front  After  Wind  has  been 

Blowing  From  Cold  to  Warm  Side  1  Hour. 
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Blowing  From  Cold  to  Warm  Side 


The  results  of  the  computations  are  given  in  Figures 
14  and  15.  When  the  wind  blows  from  the  warm  to  the  cold 
side,  AT  and  Ae  are  decreased  (Figures  10  and  11)  and  the 
corresponding  heat  exchange  components  are  reduced  (Figure 
14) .  The  opposite  is  the  case  when  wind  blows  from  the 
cold  to  the  warm  side  (Figure  15) .  The  latter  process  pro 
vides  additional  moisture  to  the  warm  side  of  the  front 
through  advection  and  the  increase  in  Qe  which  results  in 
upward  motions  and  cloud  formation. 
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HEAT  EXCHANGE  g  cal  cmw*  ( 2  4  h  ) 


DISTANCE,  km. 

Figure  14.  Sensible  and  Latent  Heat  Exchange  Across  a  Front; 
Wind  Blowing  From  Warm  to  Cold  Side. 
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HEAT  EXCHANGE  g  cal  cm'^  (2Ah) 


DISTANCE,  km. 

Figure  15.  Sensible  and  Latent  Heat  Exchange  Across  a  Front; 
Wind  Blowing  From  Cold  to  Warm  Side. 


31 


5.  NUMERICAL  SIMULATION  OF  SURFACE  PRESSURE  CHANGE, 

Vertical  motion,  and  "zusatzwind"  at  oceanic?  fronts 

v  i  -  0fO  ,,oT 

The  numerical  simulation  of  the  feedback  effects, 

briefly  reported  below,  was  aimed  at  obtaining  a  simplified 

q  font.  ,,  •!  <'.I  jvoi  ;.au,'.soiq  ixoowjou  aaenrfoirtJ  ~  ri 

semiquantitative  picture  of  the  nature  of  these  effects  at 

oxledsibu  jH  +  mot 3  oi'J-  io  erufrnqtr 

oceanic  frontal  regions. 

q  ;  pni  tnoqrjfo  ,  insdan  i  o  *  /■. 

The  equation  of  state 

P-': 


pR'T* 

P  =  - 

T 


(15) 


to  aci.^tro'. 


T 

where  T*  =  -  indicates  that  pressure  changes 

(1  -  0.379  ea/p)  *  y 

across  an  oceanic  front'  can  be34xpected  due3 to3 the  effect 
of  the  fronts  On  the  differential  heating  and/or  cooling  of 
the  lower  layers  of- the  air.  r 

The  first" law  of 3 thermodynamics  (equation  (4))  in  the 


adiabatic  form: 


dQ  =  CDdT  -  —  dp  (16) 

p  P 

indicates  about  a  2  mb.  lowering  of  pressure  during  the 
travel  of  air  from  the  warm  to  the  cold  side,  across  the 
front.  Cormier  and  Kindle  (1971) ,  using  different  ap¬ 
proaches  found  that  simulated  effects  in  their  numerical 
models  were  less  than  the  observed  effects. 

The  Sawyer  and  Bushby  (see  Sawyer,  1963)  equation  was 
used  for  approximations  of  the  vertical  velocity  (wq)  caused 
by  surface  heating: 
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(17) 


to. 


Tc„  (P, 


- Q 

Pj)  RArp 


where : 

h1  -  thickness  between  pressure  levels  P0  and  Pt 
Tp  -  departure  of  the  lapse  rate  from  the  adiabatic 
A  -  constant,  depending  on  P0  and  Pt 
R  -  gas  constant 
C  -  specific  heat 

Jr 

g  -  acceleration  of  gravity 
Q  -  heat 

T  -  mean  temperature  of  the  column. 

The  results  of  the  computations  are  shown  in  Figure  16 
and  indicate  upward  motion  when  the  wind  blows  toward  the 
warm  side  and  downward  motion  when  the  wind  is  in  the  op¬ 
posite  direction.  As  pointed  by  Sawyer  (1965)  it  is  at 
present  not  clear  how  baroclinic  disturbances  caused  by 
heat  exchange  at  the  surface  affect  the  vertically  inte¬ 
grated  flow.  The  column  of  air  is  affected  not  only  by  con¬ 
vergence  at  one  level  but  also  by  the  balancing  divergence 
at  another  level.  However,  increased  upward  motion  in  an 
area  (such  as  caused  by  air  flowing  from  the  cold  towards 
the  warm  side  of  an  oceanic  front)  will  result  in  an  in¬ 
creased  generation  of  cyclonic  vorticity  near  the  ground. 
This  may  explain  the  well-observed  fact  that  extratropical 
cyclones  tend  to  intensify  and  accelerate  near  oceanic 
frontal  regions. 
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Figure  16.  Estimate  of  Vertical  Motion  at  the  Front,  Using 

Sawyer-Bushby  Equation. 
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The  simplified  computation  of  pressure  change  across 
an  oceanic  front  cannot  be  considered  fully  applicable  to 
real  nature  but  will,  however,  give  indications  of  processes 
in  frontal  regions.  The  change  in  pressure  gradient  across 
the  front  caused  by  heating  or  cooling  of  the  surface  layer 
or  by  vertical  motion  will  cause  an  additional  wind  component 
referred  to  by  the  German  term  "zusatzwind" .  This  wind  com¬ 
ponent  was  computed  numerically  using  the  geostrophic  wind 
equation  and  the  computed  pressure  change  which  resulted 
from  heating  or  cooling  of  the  lower  layer  of  the  air.  In 
the  case  of  the  Sawyer,  Bushby  equation,  the  vertical  motion 
was  first  converted  to  pressure  change.  An  example  of  the 
zusatzwind  computed  from  the  equation  of  state  is  given  in 
Figure  17.  Areas  of  convergence  and  divergence  between  the 
wind  and  the  zusatzwind  are  shown.  If  the  wind  direction  is 
reversed,  the  area  of  acceleration  and  convergence  remain  in 
the  same  relative  positions  to  the  basic  wind  direction,  but 
change  in  relation  to  the  orientation  of  the  front  (i.e., 
the  convergence  moves  to  the  cold  water  side  of  the  front) . 
The  zusatzwind  derived  from  the  thermodynamic  equation 
(Figure  18)  is  considerably  weaker  and  the  patterns  differ. 
Which  of  the  simplified  computations  gives  the  more  realistic 
picture  is  difficult  to  decide  at  this  stage.  The  answer 
to  this  question  will  be  obtained  from  actual  observations 
and  from  the  computations  with  fine  mesh  window  type  (zoomed) 


ADDITIONAL  WIND  cm  sec 


WIND 


-2  0-1 - T - r - 

0  100 

DISTANCE,  km  . 


200 


Figure  17 . 


Zusatzwind  From  Equation 
From  Warm  to  Cold 


of  State,  Wind  Blowing 
Side. 
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ADDITIONAL  WIND  cm  sec 


Figure  18.  Zusatzwind  From  Thermodynamic  Equation;  Wind 
Blowing  From  Warm  to  Cold  Side. 


- 


models.  Attention  could  be  called  to  the  review  of  this  sub¬ 
ject  by  Roll  (1965)  who  pointed  out  the  occurrence  of  unusual 
mesoscale  winds  at  sea  which  do  not  follow  the  classical 
relationship  between  reported  surface  wind  and  the  computed 
gradient  wind. 

The  low  resolution  of  the  computational  grid  in  hemi¬ 
spheric  atmospheric  models  inhibits  the  presentation  of  meso¬ 
scale  features  such  as  oceanic  frontal  effects  and  coastal 
effects.  However,  the  sum  of  mesoscale  features  affects  the 
hemispheric  scale  feedback  and  proper  means  must  be  found  to 
include  this. 

There  are  clear  indications  in  published  data  that  the 
wind  fields  respond  to  oceanic  fronts.  One  of  the  studies 
by  Verploegh  (1954)  is  reproduced  in  Figure  19  which  shows 
a  remarkable  relationship  between  air/sea  temperature  and 
wind  speed.  Surface  winds  in  excess  of  twice  the  gradient 
wind  along  the  North  Wall  of  the  Gulf  Stream  (especially 
during  cold  air  outbreaks)  are  also  pointed  out  by  Cormier 
and  Kindle  (1971) . 

The  major  oceanic  fronts  are  usually  much  more  pronounced 
than  the  ones  assumed  in  the  numerical  example.  Considering 
the  effects  of  the  oceanic  fronts  on  the  surface  air  (as 
demonstrated  in  this  paper,  and  by  empirical  evidence  not 
fully  described  here)  it  becomes  apparent  that  the  oceanic 
fronts  are  a  significant  feature  in  the  placement  of  the  mean 
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Figure  19.  Changes  of  Meteorological  Elements  Along  the 
Track  of  M/S  "Borneo"  Off  Somali  Coast, 

August  1952. 
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positions  of  atmospheric  fronts  near  the  vicinity  of  the 
former. 

The  theories  and  results  of  numerical  computations 
should  be  verified  by  measurements  in  nature  (e.g.  frequent 
measurements  on  ships  of  opportunity  while  crossing  oceanic 
fronts) .  The  verification  of  the  results  is  important  as 
the  same  phenomena  occurs,  to  a  lesser  degree,  in  all  parts 
of  the  oceans  when  winds  blow  across  sea  surface  isotherms. 
As  the  changes  are  most  pronounced  at  oceanic  frontal  re¬ 
gions,  more  information  is  available  in  these  pronounced 
thermal  contrast  areas. 
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6.  CONCLUSIONS 


a.  The  Amot  -  Mosby  theory  of  the  response  of  surface  air 
temperature  (Ta)  and  water  vapor  pressure  (ea)  to  the 
corresponding  properties  of  the  sea  surface  is  adapted 
for  synoptic  analysis  and  prediction  of  Ta  and  ea  over 
the  oceans . 

b.  When  this  method  is  applied  numerically  to  oceanic 
frontal  region  conditions,  relatively  rapid  modifica¬ 
tion  of  surface  air  properties  through  heat  and  moisture 
exchange  results.  Several  peculiarities  of  this  modifi¬ 
cation  process  are  observed  in  relation  to  winds  and 
oceanic  fronts  (e.g.  the  "overshooting"  of  sea-air  tem¬ 
perature  difference  and  others) .  The  peculiarities  are 
physically  explainable,  however. 

c.  When  the  rate  of  change  of  sea  surface  temperature  along 
a  trajectory  of  the  air  remains  constant,  the  sea-air 
temperature  difference  reaches  an  "equilibrium  value" 

in  about  five  hours . 

d.  The  existence  of  an  "equilibrium"  difference  between  the 
surface  air  temperature  and  sea  surface  temperature  (or 
respective  water  vapor  pressures)  does  not  mean  a  no 
heat  and  moisture  exchange  condition,  but  merely  in¬ 
dicates  that  the  supply  of  heat  (and/or  moisture)  from 
the  sea  equals  its  removal  (by  turbulence)  from  the 
surface  layers  of  the  air. 
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Due  to  pronounced  differential  heating/cooling  of  the 
lower  layers  of  the  air  in  oceanic  frontal  regions , 
relatively  sharp  mesoscale  pressure  gradients  are 
created  which  will  modify  the  local  winds . 

The  mean  positions  of  atmospheric  polar  and  arctic 
fronts  over  the  oceans  are  located  near  the  corre¬ 
sponding  positions  of  the  oceanic  fronts.  This  study 
demonstrates  that  the  oceanic  fronts  are  significant 
controlling  factors  in  the  relative  mean  positions  of 
atmospheric  fronts. 

Due  to  the  sharp  horizontal  thermal  gradients  and  pro¬ 
nounced  effects  of  these  gradients  on  the  lower  atmo¬ 
sphere,  the  oceanic  frontal  regions  are  the  most  ideal 
areas  for  the  study  and  modelling  of  the  mesoscale  feed¬ 
back  from  the  ocean  to  the  atmosphere. 

As  the  mesoscale  processes  affect  the  large-scale 
features,  proper  methods  must  be  devised  for  their 
inclusion  into  the  large-scale  numerical  predictions 
of  the  atmosphere. 
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